FNAL-PUB-95/153, |h^p~-ex/9506007 



Limits on v^^) — ► v t (Pt) an d Vyj^py) —> ^ e (^e) Oscillations from a 
Precision Measurement of Neutrino-Nucleon Neutral Current 

Interactions 



K. S. McFarland 5 , D. Naples 5 , C. G. Arroyo 4 , P. Auchincloss 8 , P. de Barbaro 8 



ON 



A. 0. Bazarko 4 , R. H. Bernstein 5 , A. Bodek 8 , T. Bolton 6 , H. Budd 8 , J. Conrad 4 , 



R. B. Drucker 7 , D. A. Harris 8 , R. A. Johnson 3 , J. H. Kim 4 , B. J. King 4 , T. Kinnel 9 , 
G. Koizumi 5 , S. Koutsoliotas 4 , M. J. Lamm 5 ,W. C. Lefmann 1 , W. Marsh 5 ,C. McNulty 4 , 
S. R. Mishra 4 , P. Nienaber 10 , M. Nussbaum 3 , M. J. Oreglia 2 , L. Perera 3 , P. Z. Quintas 4 , 

> : 

A. Romosan 4 , W. K. Sakumoto 8 , B. A. Schumm 2 , F. J. Sciulli 4 , W. G. Seligman 4 , 
^ ' M. H. Shaevitz 4 , W. H. Smith 9 , P. Spentzouris 4 , R. Steiner 1 , E. G. Stern 4 , 

o : 

in ■ M. Vakili 3 ,U. K. Yang 8 , 

on ; 

>^ ' 1 Adelphi University, Garden City, NY 11530 

U : 

fS i 2 University of Chicago, Chicago, IL 60637 

<d : 

'-^ ' 3 University of Cincinnati, Cincinnati, OH 45221 

>'• 

4 Columbia University, New York, NY 10021 

5^ ■ 

5 Fermi National Accelerator Laboratory, Batavia, IL 60510 
6 Kansas State University, Manhattan, KS 66506 
7 University of Oregon, Eugene, OR 97403 

8 University of Rochester, Rochester, NY 14627 

9 University of Wisconsin, Madison, WI 53706 
10 Xavier University, Cincinnati, OH 45207 

(February 7, 2008) 



1 



Abstract 



We present limits on v^ty^) — ► v T (V T ) and f^(F^) — > ^ e (F e ) oscillations 

based on a study of inclusive vN interactions performed using the CCFR 

massive coarse grained detector in the FNAL Tevatron Quadrupole Triplet 

neutrino beam. The sensitivity to oscillations is from the difference in the 

longitudinal energy deposition pattern of u^N versus u T N or u e N charged 

current interactions. The energies ranged from 30 to 500 GeV with a 

mean of 140 GeV. The minimum and maximum za, flight lengths are 0.9 km 

and 1.4 km respectively. For —>■ v T oscillations, the lowest 90% confidence 

upper limit in sin 2a of 2.7 x 10~ 3 is obtained at Am 2 ~ 50 eV 2 . This result 

is the most stringent limit to date for 25 < Am 2 < 90 eV 2 . For — > v e 

oscillations, the lowest 90% confidence upper limit in sin 2 2a of 1.9 x 10~ 3 is 

obtained at Am 2 ~ 350 eV 2 . This result is the most stringent limit to date 

for 250 < Am 2 < 450 eV 2 , and also excludes at 90% confidence much of the 

high Am 2 region favored by the recent LSND observation. 
14.60.Pq, 13.15.+g, 12. 15. Mm 
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The mixing of non-degenerate neutrino mass eigenstates would lead to oscillations of one 
neutrino type into another. For mixing between two generations, the oscillation probability 
is 



where Am 2 is |m 2 — m\\ , a is the mixing angle, E v is the neutrino energy, and L is the distance 
the neutrino travels between production and observation. Nonzero neutrino mass and mixing 
would have important implications for cosmology and particle physics. Neutrino oscillations 
may explain observed neutrino deficits from the sun and from atmospheric sources. 



pearance through exclusive processes. For example, the FNAL-E531 limit |lj] comes from 
searching for a detached vertex from a tau decay in emulsion. A fine grained detector {i.e. 
emulsion, or a low density fine grained calorimeter such as was used by CHARM II ||) is nec- 
essary to be sensitive to low mixing angles through exclusive modes. In the case of — ► v e 
oscillations, the best limits from accelerator experiments come from fine-grained calorimetric 



(e.g.: BNL-E734 §, BNL-E776 @) or fully active detectors (e.g.: KARMEN §, LSND [§, 



searching for quasi-elastic charged-current production of electrons. The LSND experiment, 
using a liquid scintillator neutrino target, has recently reported a signal consistent with 
—>■ Pg at a sin 2 2a of ~ 1CT 2 and a Am 2 ~ 1 eV 2 @. 

In this report, we obtain results with comparable sensitivity to E531, CHARM II for 
V/j, — > v T oscillations and to BNL-E734, BNL-E766 and KARMEN for — > v e oscillations 
at Am 2 ~ 40 eV 2 using the massive and relatively coarse grained CCFR detector. The 
main advantage of this type of detector is increased interaction probability which will be 
particularly important in a low flux, long baseline neutrino beam |7| . Our result establishes 
the sensitivity of such detectors to small mixing angles. 

The CCFR detector || consists of an 18 m long, 690 ton target calorimeter with a mean 
density of 4.2 g/cm 3 , followed by an iron toroid spectrometer. The target calorimeter consists 
of 168 iron plates, 3m x 3m x 5.1cm each. The active elements are liquid scintillation 




(1) 



To date the best limits for v. 



v T oscillations are derived from searches for v T ap- 
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counters spaced every two plates and drift chambers spaced every four plates. There are a 
total of 84 scintillation counters and 42 drift chambers in the target. The toroid spectrometer 
is not directly used in this analysis. 

The Tevatron Quadrupole Triplet neutrino beam is created by decays of pions and kaons 
produced when 800 GeV protons hit a production target. A wide band of secondary energies 
is accepted by focusing magnets. The production target is located about 1.4 km upstream 
of the neutrino detector. The production target and focusing train are followed by a 0.5 km 
decay region. The resulting neutrino energy spectra for u^, F M , u e , and z/ e at the detector are 
shown in Figure p]. The beam contains a 2.3% fraction of electron neutrinos and a negligible 
fraction of tau neutrinos (less than 10 -5 ) which result primarily from D s decay. 

Neutrinos are observed in the target calorimeter via their neutral current and charged 
current interactions, charged current events are characterized by the presence of a muon in 
the final state which deposits energy in a large number of consecutive scintillation counters 
as it travels through the calorimeter. Neutral current events have no muon and deposit 
energy over a range of counters typical of a hadronic shower (5 to 20 counters) . Accordingly, 
we define "short" events as those which deposit energy over an interval of 30 or fewer 
scintillation counters. The ratio R30 is defined to be the number of short events divided by 
the number of long events ||. This ratio is strongly dependent on the ratio of neutral to 
charged current events which is a function of the electroweak mixing angle, sin 2 6w 

Assuming the validity of the Standard Model, sin 2 6\y is accurately measured from other 
processes. We can use these measurements to predict the ratio of neutral to charged current 
events in the CCFR detector, and thus R30. The presence of v T or additional v e in the 
neutrino beam would cause the measured R30 to be larger than its calculated value because 
most charged current tau and electron neutrino interactions do not produce a muon in the 
final state and will thus appear "short" . 

In this study, we attribute any deviation in our measured R30 from the predicted value 
to — > v T or — > v e oscillations. This technique, which has been discussed previously 
in the literature |7j] [[!(]] [|TTJ, assumes that only one of the two types of flavor oscillation 
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FIG. 1. Neutrino energy spectra for u^, u^, v e , and V e at the CCFR detector for the FNAL 
wideband neutrino beam. (Monte Carlo based on measured relative and fluxes). 
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contributions to a change in R^o, and is therefore conservative since both types of oscillation 
would increase the measured R 30 . 

We used a detailed Monte Carlo to relate a given — ► v T or — > i/ e oscillation 
probability to the quantity H30. The sin 2 valueQ in the on-shell renormalization scheme 
of 0.2232 ± 0.0018 is input to the Monte Carlo. The other inputs to the Monte Carlo are 
parameterizations of the measured CCFR detector responses |§, nucleon structure functions 
|16| , and relative neutrino beam fluxes extracted from the charged current data sample [jlT 



The v e flux is modeled in a detailed beamline simulation, normalized by the observed 
flux |J. The same beamline simulation is used to tag the decay location for each pion and 
kaon and thus the creation point of each along the beamline. The measured flux gives 
the number of z/^'s at the detector. Piv^ — > v T ,e) is determined from Eq. [I] and the beamline 
simulation. We assume P{v^ — > v T)e ) = P(?n ~ > &V,e) ( a consequence of CP invariance). 
The number of v^s produced in the beamline is then the number observed at the detector 
divided by 1 — Piy^ — > z/ r e ). The predicted electron neutrino flux is rescaled to the produced 
number of i^'s. The tau or electron neutrino flux from neutrino oscillations is calculated by 
multiplying the produced number of v^s by Piv^ — > v T) e)- 

To simulate v T ^ e interactions in our detector we assumed the v Tfi neutral current cross 
section is the same as for interactions. The v T charged current cross section was calcu- 



lated including mass suppression terms. Following [18] we used the approximation that the 



structure functions _F 4 = 0, and xF 5 = 2xFi. The kinematic suppression for massive particle 
production was also taken into account. The Monte Carlo program TAUOLA [|Tj| was used 



1 This value for sin 2 9w is obtained using the world average value My/ measurement [12|, the 
prediction from the measured Mz, and the average of all LEP and SLD Z-pole measurements from 



13fl . The Mz extraction is corrected for the recent re-evaluation of o.em{M z ) by Swartz 14]. A 



top mass of 180 ± 12 GeV [15| and 60 < Mgiggs < 1000 GeV are used to convert from the MS and 



Mz schemes to the on-shell scheme used here. 
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FIG. 2. i?3o as a function of E ca i for the data (points). The filled band shows Monte Carlo 
assuming no oscillations with la systematic errors added in quadrature. Data points show sta- 
tistical errors only. The dotted curve corresponds to Vn — > v T oscillations with Am 2 = 3500 eV 2 
and sin 2 2a = 0.05 and the dashed curve to Am 2 = 100 eV 2 and sin 2 2a = 0.05. The curve for 
Am 2 = 100 eV 2 corresponds to lower energy neutrinos for which less energy is deposited in the 
calorimeter. At high Am 2 the high E ca i events are most sensitive to v T or v e appearance. 

to simulate tau decays. We define E ca i as the energy deposited in the calorimeter in the first 
twenty counters following the event vertex. For v T>e charged current events E ca i includes 
the visible energy from the tau decay. Events are required to deposit a minimum energy 
of 30 GeV in the target calorimeter. The contributions from quasi-elastic and resonance 
production are suppressed by this requirement. 

Events were selected using a calorimeter trigger fully sensitive for E ca \ above 20 GeV. 
To ensure event containment, the fiducial volume of the detector is limited to a central 
cylindrical region 30" in radius and excludes events which began in the first 6 counters or 
the last 34 counters of the detector. The resulting data sample consisted of about 450,000 



Source of Error 



Am 2 = 3500 eV 2 



310 eV 2 



80 eV 2 



statistical 

v e beam content 

detector systematics 

charm mass 

charm sea 

sin 2 9w 

other model 

Total 



2.4 x 10~ 3 

2.5 x 10" 3 
2.2 x 10" 3 
2.4 x 10~ 3 
1.2 x 10~ 3 
1.7 x 10" 3 
0.6 x 10~ 3 
5.2 x 10~ 3 



1.8 x 10~ 3 
2.0 x 10~ 3 

2.0 x 10~ 3 
1.6 x 10~ 3 
0.8 x 10~ 3 
1.2 x 10~ 3 
0.6 x 10~ 3 

4.1 x 10~ 3 



2.1 x 10~ 3 
1.4 x 10~ 3 
2.0 x 10~ 3 

2.4 x 10~ 3 
1.0 x 10~ 3 

1.5 x 10 -3 
0.7 x 10~ 3 
4.4 x 10~ 3 



TABLE I. The change in sin 2 2a for z/„ — > v T from a one sigma shift in the dominant errors. 
The row labeled "total" includes these and many smaller uncertainties added in quadrature. 

events. The data and Monte Carlo are divided into 21 E ca i bins. For each Am 2 , the Monte 
Carlo prediction for R 3Q (E ca i, sin 2 2a) is compared with R 30 (E ca i) from the data. Figure 
H] shows the R30 distribution as a function of E ca \ for the data and for the Monte Carlo 
simulation. The detailed shape of R 30 (E ca i) depends on many competing effects which are 
put into the Monte Carlo, but is dominated by the variation of short charged current events 
with E ca i and by the contribution from the predicted u e flux. 

There are four major uncertainties in the comparison of i?3o(-E ca i) from the Monte Carlo 
to the data: the statistical error in the data, the uncertainty in the effective charm quark 
mass for charged current charm production, the uncertainty in the incident flux of u^s on the 
detector, and the uncertainty in the on-shell weak mixing angle from outside measurements. 
Other sources of systematic uncertainty were also investigated |§. Table | shows the effect 
of the uncertainties for three choices of Am 2 . 

The charm mass error comes from the uncertainty in modeling the turn-on of the charm 
quark production cross section. The Monte Carlo uses a slow-rescaling model with the 
parameters extracted using events with two oppositely charged muons in this experiment 
pOfl . This error dominates the calculation of R 3 o at low E v (and low E ca {) where the threshold 



S 



suppression is greatest. The v e flux uncertainty has a large effect on i? 30 because almost 
all charged current v e events are short events. Therefore, the relatively small (4.2% |5|) 
fractional uncertainty in the v e flux is a large effect, particularly at high E ca \ since most v e 
charged current interactions deposit the full incident neutrino energy into the calorimeter. 
This 4.2% is dominated by a 20% production uncertainty in the Kl content of the secondary 
beam which produces 16% of the v e flux. The bulk of the v e flux comes from decays, 
which are well-constrained by the observed spectrum from K^ 2 decays |§. 

The data are fit by forming a x 2 which incorporates the Monte Carlo generated effect of 
oscillations, and statistical and systematic uncertainties. A best fit sin 2 2a is determined for 
each Am 2 by minimizing the x 2 as a function of sin 2 2a and the 33 systematic coefficients, 
Cj. Best fit values of sin 2 2a with one sigma errors from the fit are shown in Tables [II] and [II. 



At all Am 2 , the data are consistent with no observed z/ M — > v r e oscillation. The statistical 
significance of the best-fit oscillation at any Am 2 is at most 1.2 sigma. 

The frequentist approach |21| is used to set a 90% confidence upper limit for each Am 2 . 



The limit in sin 2 2a at each Am 2 corresponds to a shift of 1.64 units of x 2 (sin 2 2a) from 
the minimum x 2 (at the best fit value in Table [TT]). The x 2 value for the no-oscillations case 
is 15.7/21 dof. The 90% confidence upper limit is plotted in Figure |3| for z/ M — > v T . The 
best limit of sin 2 2a < 2.7 x 10~ 3 is at Am 2 = 50 eV 2 . For sin 2 2a = 1, Am 2 > 1.4 eV 2 
is excluded, and for Am 2 > 1000 eV 2 , sin 2 2a > 8.1 x 10~ 3 is excluded at 90% confidence. 
For — > v e oscillations, the 90% confidence upper limit is also shown in Figure [3[ The 
best limit of sin 2 2a < 1.9 x 10~ 3 is at Am 2 = 350 eV 2 . For sin 2 2a = 1, Am 2 > 1.6 eV 2 is 
excluded, and for Am 2 > 1000 eV 2 , sin 2 2a > 3.8 x 10" 3 

This result demonstrates sensitivity to low mixing angles in a high mass, coarse grained 
sampling calorimeter and has implications for proposed long-baseline experiments . How- 
ever, a detailed Monte Carlo study of the sensitivity of those experiments must be performed 
to correctly apply this result. The lower energy and the lower level of statistics in the long 
baseline experiments will result in less statistical sensitivity, while having both a near and 
far detector would reduce many of the other sources of uncertainty [|ll|] listed in Table |. 
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Am 2 (eV 2 ) 


Best Fit 


Sigma 


Am 2 (eV 2 ) 


Best Fit 


Sigma 


2.0 


-2.1114 


2.0192 


185.0 


0.0050 


0.0042 


3.5 


-0.6982 


0.6676 


200.0 


0.0047 


0.0042 


5.0 


-0.3419 


0.3268 


220.0 


0.0040 


0.0041 


6.0 


-0.2373 


0.2267 


240.0 


0.0033 


0.0041 


8.0 


-0.1351 


0.1296 


275.0 


0.0022 


0.0040 


10.0 


-0.0872 


0.0838 


295.0 


0.0018 


0.0040 


15.0 


-0.0397 


0.0385 


310.0 


0.0016 


0.0040 


20.0 


-0.0229 


0.0224 


350.0 


0.0010 


0.0041 


35.0 


-0.0084 


0.0090 


400.0 


0.0002 


0.0043 


42.0 


-0.0061 


0.0069 


430.0 


-0.0002 


0.0044 


50.0 


-0.0045 


0.0056 


500.0 


-0.0002 


0.0047 


60.0 


-0.0031 


0.0047 


550.0 


0.0004 


0.0048 


70.0 


-0.0021 


0.0041 


600.0 


0.0013 


0.0050 


80.0 


-0.0013 


0.0038 


650.0 


0.0020 


0.0051 


90.0 


-0.0006 


0.0036 


700.0 


0.0026 


0.0051 


100.0 


0.0001 


0.0035 


750.0 


0.0028 


0.0050 


110.0 


0.0008 


0.0035 


800.0 


0.0027 


0.0049 


120.0 


0.0015 


0.0036 


1000.0 


0.0017 


0.0049 


135.0 


0.0027 


0.0037 


2000.0 


0.0018 


0.0050 


150.0 


0.0038 


0.0039 


3500.0 


0.0018 


0.0049 


175.0 


0.0049 


0.0041 


10000.0 


0.0018 


0.0049 



TABLE II. The result for sin 2 2a from the fit at each Am 2 for f„ — > v T oscillations 
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Am 2 (eV 2 ) 


Best Fit 


Sigma 


Am 2 (eV 2 ) 


Best Fit 


Sigma 


2.0 


-0.7856 


1.0714 


200.0 


0.0008 


0.0021 


3.0 


-0.3552 


0.4825 


225.0 


0.0001 


0.0020 


4.0 


-0.1993 


0.2720 


250.0 


-0.0001 


0.0019 


5.0 


-0.1275 


0.1746 


275.0 


-0.0001 


0.0017 


7.0 


-0.0650 


0.0897 


300.0 


-0.0001 


0.0016 


9.0 


-0.0392 


0.0548 


350.0 


-0.0001 


0.0015 


10.0 


-0.0317 


0.0448 


400.0 


0.0000 


0.0016 


20.0 


-0.0074 


0.0122 


450.0 


0.0004 


0.0018 


30.0 


-0.0027 


0.0061 


500.0 


0.0009 


0.0020 


40.0 


-0.0009 


0.0039 


600.0 


0.0015 


0.0024 


50.0 


-0.0001 


0.0028 


700.0 


0.0014 


0.0025 


60.0 


0.0004 


0.0023 


800.0 


0.0007 


0.0023 


70.0 


0.0007 


0.0019 


1000.0 


0.0009 


0.0023 


80.0 


0.0009 


0.0018 


1500.0 


0.0009 


0.0022 


90.0 


0.0012 


0.0017 


2000.0 


0.0009 


0.0023 


100.0 


0.0014 


0.0017 


5000.0 


0.0009 


0.0023 


125.0 


0.0022 


0.0018 


10000.0 


0.0009 


0.0023 


150.0 


0.0024 


0.0020 


20000.0 


0.0009 


0.0023 


175.0 


0.0017 


0.0021 


50000.0 


0.0009 


0.0023 



TABLE III. The result for sin 2 2a from the fit at each Am 2 for Vn — > v e oscillations 
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10 3 10 2 sin 2 2a 10 1 1 

FIG. 3. Excluded region of sin 2 la and Am 2 for f„ — > v T ,e oscillations from this analysis at 
90% confidence is shown as dark, solid curves. 



In conclusion, we have used a new analysis method to search for — > z/ T ,e oscillations 
with a coarse-grained calorimetric detector. We see a result consistent with no neutrino 
oscillations and find 90% confidence level excluded regions in sin 2 2a- Am 2 space. This 
result is the most stringent limit to date for — > v T oscillations with 25 < Am 2 < 90 eV 2 
and for v„ — >■ v e oscillations with 250 < Am 2 < 400 eV 2 . 
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